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ABSTRACT 

Context. Stray light caused by scattering on optical surfaces and in the Earth's atmosphere degrades the spatial resolution of ob- 
servations. Whereas post-facto reconstruction techniques are common for 2D imaging and spectroscopy, similar options for slit- 
spectrograph data are rarely applied. 

Aims. We study the contribution of stray light to the two channels of the POlarimetric LIttrow Spectrograph (POLIS) at 396 nm and 
630 nm as an example of a slit-spectrograph instrument. We test the performance of different methods of stray-light correction and 
spatial deconvolution to improve the spatial resolution post-facto. 

Methods. We model the stray light as having two components: a spectrally dispersed component and a "parasitic" component of 
spectrally undispersed light caused by scattering inside the spectrograph. We used several measurements to estimate the two contribu- 
tions: a) observations with a (partly) blocked field of view (FOV), b) a convolution of the FTS spectral atlas, c) imaging of the spider 
mounting in the pupil plane, d) umbral profiles, and e) spurious polarization signal in telluric spectral lines. The measurements with a 
partly blocked FOV in the focal plane allowed us to estimate the spatial point spread function (PSF) of POLIS and the main spectro- 
graph of the German Vacuum Tower Telescope (VTT). We then used the obtained PSF for a deconvolution of both spectroscopic and 
spectropolarimetric data and investigated the effect on the spectra. 

Results. The parasitic contribution can be directly and accurately determined for POLIS, amounting to about 5 % (0.3 %) of the (con- 
tinuum) intensity at 396 nm (630 nm). The spectrally dispersed stray light is less accessible because of its many contributing sources. 
We estimate a lower limit of about 10 % across the full FOV for the dispersed stray light from umbral profiles. In quiet Sun regions, 
the stray-light level from the close surroundings (d < 1") of a given spatial point is about 20 %. The stray light reduces to below 
2% at a distance of 20'' from a lit area for both POLIS and the main spectrograph. The spatial deconvolution using the PSF obtained 
improves the spatial resolution and increases the contrast, with a minor amplification of noise. 

Conclusions. A two-component model of the stray-light contributions seems to be sufficient for a basic correction of observed spectra. 
The instrumental PSF obtained can be used to model the off-limb stray light, to determine the stray-light contamination accurately for 
observation targets with large spatial intensity gradients such as sunspots, and also to improve the spatial resolution of observations 
post-facto. 

Key words. Sun: chromosphere - photosphere - Methods: data analysis - Line: profiles 



1. Introduction 



The importance of stray-light contributions to observed spec 
/\ tra and images was realized ve ry early on and continues to 
^ [ be a problem today (see, e.g.. iHenouxl Il969l: iMattid 1971 



Martinez Pillet''l992^:"Chaeetal''l998':" Wedemeve^ 
DeForest et al. 2009; Mathew et al. 2009|). The intrinsically lim- 
ited optical quality of the reflective surfaces of mirrors or of the 
glass of lenses leads to scattering of photons in the light path. 
In addition, scattering during the passage through the Earth's at- 
mosphere by, for instance, dust particles spreads the light from 
the full solar disk to any point of an observed restricted field of 
view (FOV). Depending on the wavelength and the target of the 
observations, the stray light can amount to a significant fraction 
of the observed intensity and can have a strong impact on the fi- 
nal results of the data analysis. Therefore, many inversion codes 
that determine ph ysical quantities from observed spec tra, such 
as the SIR code (iRuiz Cobo & del Toro Iniestalll992h . include 
an explicit treatment of this contamination by using a separately 



provid ed stray-light profile (see also, e.g.. lOrozco Suarez et al.l 
12007^ . 

For analysis techniques that make direct use of the ob- 
served spectra or of broad-band images, the stray light has to 
be dealt with in advance. This problem had to be tackled in 
the context of the accurate determination of sunspot intensi- 
ties to derive their temperature at continuum forming layers 
(Kneer&Mattig 1968; Maltbv & Mvkland|ll969|; iMaltbvll 19701: 
iMattig..l971.:.Tritschler & Schmidt200Z) or in the determination 
of the continuum contrast of solar granulation ([Mathew et"aD 
120091: IWedemever-Bohm & Rouppe van der VoortI l2009l) . It is 



also important in studies of solar chemical abundances that 
seek to derive consistent results from d ifl^erent spectral lines 
(lAsplund et al.ll2009l:lF^ian et al.ll2010l) . 

The stray-light contribution is often decomposed in difl'er- 
ent sources: 1. the atmospheric stray light caused by large-scale 
scattering in the Earth's atmosphere, with a slow temporal evo- 
lution because of the change of air-mass in the light path during 
the day; 2. the "blurring", i.e., the rapidly changing amount of 
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stray light caused by the fluctuations of the refractive index in 
the Earth's atmosphere (which are nowadays commonly referred 
to as "seeing"); and 3. instrumental stray light, i.e., scattering 
in the telescope/instrument optics or (static) stray-light efl'ects 
of the telescope/instrument caused by the geometry of the fi- 
nite aperture and t he possible presence of ( central) obscurations 
in the li ^ht path (Zwaan "1965; "Staveland'" 19701; iMattigl 119831: 
iMartinez Pillet 1992; Wedeme y er-B ohm 2008). For space-based 
observations, naturally, the first two points are absent. 

The stray-light contamination is quantified in the "spread 
function", whose large-scale and small-scale contributions are 
usually approximated as analytical functions of varying shapes, 
e.g., Gau ssian or Lorentzian, or with combinations of a number 
of them (lMattidll97ll: IWedemever-Bohmll2008l: [PeForest et aP 



servations, whether broad-band imaging or narrow-band spec- 
troscopy or spectropolarimetry with suitable 2D instruments 
(see, e.g., Tritschler et al. 2002; Mikurda et al. 200 6; Cavallinj 
20061; IPuschmann et al.ll2006l: IScharmer et aL[l2008nBeck et al 



l2009l:lMathew et al.ll2009h 7 For describing the purely instrumen- 
tal eff'ects without temporal dependence and in night-time as- 
tronomy, the label "spatial point spread function" (PSF) is com- 
monly used because the function describes the spatial shape that 
the light from a point source (as a star can be approximated to 
be) would attain after passing through the optical system. 

To determine the exact shape of the spread function, ei- 
ther theoretical calculations of the optical systems or observa- 
tions with an only partly illuminated FOV are used. The theo- 
retical calculations are usually limited to the time-invariant ef- 
fects of the optics, while for the varying contributions from, 
e.g., the seeing only some statistically averaged efl'ects may 
be assumed. For determining the spread function in the so- 
lar case, suitable observations can be obtained in difl'erent 
ways jM attig 197H lBriand et al.ll2006l:IWedemever-B'B hm 2008; 
iDeForest et al.l 120091) . namely by planetary transits in front of 
the Sun, by observations near and ofl' the solar limb, or by ar- 
tificial occulters in the light path. Planetary transits are actually 
the best suited because they provide the perfect reference of an 
intrinsically zero light level, but the diameters of the planets are 
usually smaller than the extent of the PSF and thus do not al- 
low to determine the far wings of the PSF; moreover, transits are 
rare events. Thus, the radial variation of the residual intensity 
beyond the solar limb ("aureole") has often been used to deter- 
mine the spread function of ground-based telescopes. Some of 
the older observations, however, have to be treated with some 
care because they were done with telescopes without an aper- 
ture field stop; i.e., the telescopes created a full-disk solar im- 
age in the focal plane without a restriction to a limited FOV. 
This implies that there was a large contribution from the full- 
disk stray light in those observations. Additionally, these data 
were usually obtained without active compensation for seeing 
eff'ects by a corr elation tracker (iBalle steros et al. 19 9^ or adap- 
tive optics (AO , Ivon der Liihe et al.ll2003l: IScharmer^ al.ll2003l: 
lRimmeldl2004l), there fore mixing seeing and atmospheric scat- 
tering. iMattigl (1 19831) found that the stray light from near the 
solar limb (up to one solar radius above it) is dominated by in- 
strumental and seeing efl'ects. The atmospheric scattering only 
becomes dominant at larger distances to the Sun or in the case 
of coronagraphs, where the illumina tion from the solar disk is 
blocked inside the telescope optics. iBriand et"an (l2006l) show 
that the determination of the spread function from the aureole 
does not necessarily provide a good measure of the instrumen- 
tal stray light because the telescope/instrument optics are only 
partly illuminated near the limb, which yields a difl'erent amount 
of scattering than for the full illumination on disk center. 

Even if the static instrumental part of the PSF is known 
from either calculations or measurements, the time- variant part 
caused by the Earth's atmosphere is usually unknown in the 
case of ground-based observations. For intrinsically 2D ob- 



820101) , the instantaneous PSF can be estimated using series of 
rapid short-e xposed images. Two difl'erent approaches, speckle 
imaging (e.g., von der Liihel 19931) and blind deconvolution (e.g., 
t^an Noort et al. 2005, and references therein), then allow a post- 
facto reconstruction of the most probable original solar image 
that lies behind the observed image series. The reconstruction 
algorithms are, however, only capable of recovering the instan- 
taneous PSF within some limi ts and are also p artially blind to 
the static parts of the PSF IScharmer et al.ll201Ql) . 

For slit- spectrograph observations, at any given moment of 
time, only a ID slice of the solar surface is available, whereas 
the PSF is a 2D function that therefore also introduces stray 
light from locations that are not covered by the slit in that mo- 
ment. Thus, even if both the static and the time- variant parts of 
the PSF are known, it is not straightforward to use this informa- 
tion to im prove th e spatial resolution of slit- spectrograph data 
post- fa cto. iKeller & Johannesson ( 19^5b andlSiitterlin & Wieh^ 
(120001) introduced a method similar to a speckle deconvolution, 
requiring a fast series of slit-spectra. In our contribution, we 
tested the application of a post-facto correction of "standard" 
slit- spectrograph data for the static part of the instrumental PSF 
as obtained from a dedicated set of measurements. 

The stray-light contribution becomes significant when the 
light level in observations is low, as for observations of sunspots, 
near the solar limb, or for spectroscopy of very deep spec- 
tral lines such as Ha or Can H. For a recent study of the 
center-to-limb variation (CLV) of Ca ii H spe ctra taken with the 
POlarimetric LIttrow Spectrograph (POLIS; iBeck et al.l l2005h 
at the German Vacuum Tower Telescope (VTT, ISchroter et al.l 
Il985l) . we thus had to investigate the stray-light contamination 
in detail. During the CLV run and some of the subsequent obser- 
vational campaigns, we obtained various data sets for character- 
izing the stray-light contamination for both POLIS and the main 
spectrograph of the VTT. Because the measurements also allow 
one to estimate the spatial PSF of the instrument, we addition- 
ally investigated the option of using the known PSF for a spatial 
deconvolution of the slit- spectrograph observations. 

In Sect. [21 we outline the basic approaches to correcting 
stray light based on the observed spectra and generic estimates 
at first (Sect. 12.11) , and then explicitly taking the instrumental 
PSF into account (Sect. 12.21) . The measurements used to derive 
the generic stray-light estimates and the PSF are described in 
Sect.O Section [4] shows examples of the application of both the 
stray-light correction and the deconvolution to a few data sets. 
The application to observational data allowed us to cross-check 
and improve the stray-light estimates. Section [5] shows the rela- 
tion between the generic stray-light correction and the explicit 
use of the PSF for a few data sets. Our findings are discussed in 
Sect.[6j Section [7] presents the conclusions. 

2. Theoretical description of tlie stray-light problem 

The transfer through the optical system composed of the (fluc- 
tuating) atmosphere of the Earth, the telescope, and finally the 
instrument modifies the real intensity distribution /tme and yields 
instead an observed intensity /obs that is related to /tme by 



/obs(0 = PSF(0(8)/tme(0, 



(1) 



where PSF(t) is the instantaneous optical point spread function 
and denotes a convolution. 
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POLIS slit to camera, top view 

Fig. 1. Sketch of the light path at the VTT (drawing not to scale). 
Focal planes are denoted by blue rectangles, pupil planes by 
brown ones. Red rectangles denote an artificial blocking of the 
light path. The 45-deg mirror behind the relay optics was used to 
direct the light to either the PCO setup or towards POLIS. The 
design of POLIS is shown in side view up to the location of the 
slit (upper panel at lower right) and in top view for the rest of 
the optics (lower panel at lower right). 



The three contributions to the PSF are explained in the intro- 
duction above. For observations with real-time correction by an 
AO system, the atmospheric part is already (partially) corrected 
for, leaving the large-scale atmospheric scattering and the static 
contributions from optics. The PSF spills the light from one loca- 
tion into its surroundings and vice versa, which adds additional 
intensity to all points in the FOV via false light originating in 
other locations in the FOV, the so-called stray-light contamina- 
tion. To apply a correction for the stray-light contamination, we 
first derive a direct approach based on only the observed spectra 
and a generic estimate of the stray-light fraction in a simplified 
approach in the next section, and a more explicit model taking 
the instrumental PSF into account in Sect. 12. 2[ 

2.1. Modeling of stray-light contributions 

For the direct correction of the stray-light contamination, we 
model the intensity spectrum /obs(^, ^^y) on a certain CCD row 
in the focal plane of the POLIS CCD cameras (focal plane F4, 
see Fig.[T]) as 

/obs(/^, X,y) = /true(/^, X,y) -F Qfi Ig\oha\(^) + (^2 hoca\(^) + C0nst.,(2) 

with 

const. = dc -F parasitic light = dc -F 0^3 < /locai > ^ f(^) , (3) 



where "dc" denotes the dark current of the CCD. /tme(^, x,y) is 
the spectrum that emerges from the spatial location (x, y) on the 
Sun that corresponds to the CCD pixel row y when the slit is 
located at position x on the solar image. The average intensity 
spectrum before the first focal plane Fl is denoted by /global (^); 
the full solar disk contributes to it. Because of the CLV of the 
solar intensity, /global (^) will be dominated by contributions from 
the disk center and its surroundings, and its shape is independent 
of the telescope pointing. This term corresponds to the integra- 
tion over th e full solar dis k weighted by a spread function, as 
done in, e.s., lMattid(ll971h . It can be assumed that the stray-light 
level proportional to /global (^) is uniform across the FOV, with no 
direct relation to the solar location (x, y) where /tme(^, y) orig- 
inates. Scattering by any of the optical elements in front of Fl 
will only need to produce slight beam deviations from the optical 
axis to spill light from the full aperture to each location on the 
focal plane because of the long path length before Fl is reached. 
The factor ai then refers to all nine optical elements in front of 
Fl (7 mirrors and 2 glass windows, see the upper part of Fig. [T]). 

The average intensity spectrum behind Fl, /locai(^), comes 
only from a restricted FOV and not from the full solar disk be- 
cause of the field stop in the filter wheel near Fl. The shape of 
/local (^) depends on the telescope pointing. For the CLV obser- 
vations, we used an aperture field stop with a diameter of 40 mm 
(^^ISO'O- The diameter of the field stop in F2 inside of POLIS is 
similar to the one of the field stop in Fl. The stray light propor- 
tional to /local (^) actually is made up of two contributions: 



a^2/local(^) = 0:4 



-F < PSFinstr(Ax, Ay) X I(A, Ax, Ay) >(Ax,Ay) ^ 



(4) 



where A is the area of the field stop in Fl, and (Ax, Ay) the dis- 
tance between (x, y) and all other points inside the FOV that pass 
through the field stop. 

The first term is caused by the fact that the optical elements 
between Fl and F4 will uniformly mix the light across the full 
aperture of the field stop (as was the case for the optics upfront of 
Fl), especially because there are three additional pupil planes in- 
between. The second term depends on the spatial PSF of POLIS 
that introduces a strongly field-dependent stray-light contribu- 
tion to the spectrum from locations close to (x, y) when (Ax, Ay) 
is smaller than about 10'' (see Sect. 13.41) . PSFinstr only refers to 
the optics behind Fl. We assume that both contributions on the 
right-hand side of Eq. (|4]) can be modeled as being proportional 
to the same profile /locai(^) because in quiet Sun regions averag- 
ing over a small (a few arcsecs^) or large (several ten arcsecs^) 
area inside the local FOV will yield a similar profile. The fac- 
tor Qf2 then refers to all optics between Fl and F4 (for POLIS: 15 
mirrors, modulator, grating, interference filter, and the polarizing 
beam splitter for the 630 nm channel). The stray-light contribu- 
tions proportional to /global (^) and /locai(^) are both assumed to 
pass through a grating, and thus to be spectrally dispersed. 

The "stray light" without wavelength dependence, the con- 
stant term in Eq. (O, again contains two contributions: the usual 
dark current (dc) of the CCD cameras and a second contribution 
that corresponds to light that reaches the CCD without passing 
through the grating, called "parasitic" light. The source of the 
parasitic light photons in the case of POLIS is indicated schemat- 
ically in Fig.O which shows the optical design of POLIS just in 
front of the CCD cameras. In the light path behind the grating, 
a small pick-up mirror is used to deflect the blue part (indicated 
by a blue ray in the figure) of the dispersed spectrum towards 
the vertically mounted Call H CCD camera ("blue channel"). 



3 



C. Beck, R. Rezaei & D. Fabbian: Stray-light contamination and spatial deconvolution of spectrograph observations 



Light 
blocked 



Camera 
mirror 



Scattered 
light 




Dispersed 
spectrum 



Call H 
F4 

IF 

Pickup mirror 

3|fa 

if' 630 nm 



Fig. 2. Design of POLIS close to the CCD cameras in a side 
view (not to scale). "IF": order selecting interference filters. The 
thin red rectangle indicates the location where the light path was 
blocked for measuring the parasitic stray light. 

The part of the spectrum at a longer wavelength (indicated in 
the figure by a red ray) instead passes below the pick-up mirror 
towards the channel at 630 nm ("red channel"). Spectrally undis- 
persed scattered light inside the instrument (indicated in the fig- 
ure by black arrows) can reach the 630 nm channel CCD only 
from a small solid angle, whereas most of the scattered light that 
hits the pick-up mirror can still reach the Ca II H CCD at some 
place. Because the solar light level in the near-UV is rather low, 
especially in the line core of Ca II H, additional photons can eas- 
ily lead to a detector count rate that is significant with respect to 
the true solar spectrum. 

The parasitic stray light is generated inside of POLIS. Its 
light level should therefore scale with the amount of light that 
enters through the field stop in F2, which can be related to the 
average of < /locai > • The coefficient ^3 in Eq. ([3]) does not refer 
to any specific optical element. 

Equation Q can be simplified for any practical applica- 
tion of a stray-light correction. The full aperture contributing 
to /locai(^) is rather large, covering more than one supergran- 
ular cell and therefore providing a large-scale average profile 
similar to a full-disk spectrum, and there are no means to de- 
termine /global (^) from the observed spectra. We thus substituted 
4iobai(^) with /locai(^), ^vcn if for observations off' the disk center 
there should be some difference between the two spectra because 
/local (^) changes with the telescope pointing. To apply the correc- 
tion to observations, we here implicitly assume that the average 
profile /local (^) is always calculated from quiet Sun regions with 
a granular pattern. In the case of, e.g., sunspot observations, this 
implies that only a subfield of the observed FOV is used for the 
average. Because the dark current can be directly measured and 
subtracted, the simplified version of Eq. ^ then reads as 



/obs(/^, X,y) = /true(^, X,y) + a/iocal(^) < /local > 



(5) 



Even if the equation depends solely on two parameters, it turns 
out that determining a is far from straightforward. In Sect.O we 
describe the results of our several different approaches to deter- 
mine a and JS. 



2.2.1. First-order correction 

To derive a first-order correction for the stray light caused by the 
instrumental PSF, we reformulate and discretize Eq. ([T]) as 



hhs(x,y) = /tme(-^, j) + ^ K(X -X\y- y)Iimc(x\y) 

x',y' 

-Y,K(x-x\y-y)h^,(x,y) (6) 

x',y' 

= (1 - a)Iiruc(x. J) + ^ /^(X - J - y )/true(x\ y ) , (7) 

x',y' 

where K is the kernel describing the instrumental PSF and the 
sums are to be executed for all x\y' ^ x,y. 

The last two terms on the right-hand side of Eq. © describe 
the stray light introduced into a location (x,y) from other lo- 
cations in the FOV ("gain" term) and the "loss" of intensity to 
them, respectively. It can be shown that the first-order correction 
for the gain term is given by 

/tme(-^, y) = /obs(-^, ~ X -^'^y- y )/obs(-^', y ) , 



(8) 



x',y' 



when all terms proportional to a, a^,Ka, are neglected (see 
Appendix lAl). Equation ^ thus corresponds to the approach of 
calculating the stray light that enters at (x,y) as the observed 
intensity in the surroundings weighted by a PSF with a form as 
in the second term on the right-hand side of Eq. (|4]). 

For an application of the correction, one has to take into ac- 
count that for slit- spectrograph observations the scanning direc- 
tion (denoted by x) and the direction along the slit (denoted by 
y) are not fully equivalent. All spectra lohsi^^y') for x' x are 
taken at a different time t' , and only the profiles in the individ- 
ual spectrum /obs(-^,j') are obtained simultaneously for all y' . 
Because the solar scene is changing with time during the scan- 
ning, it can therefore be necessary to restrict the sum in x' to the 
temporally "close" surroundings instead of using the full range 
of the known PSF. The time lag between two scan steps is given 
by the integration time for each step ^mteg times their distance Ax, 
where Ax • ^integ should be significantly smaller than the typical 
granular time scale of about five minutes. 

The loss term in Eq. d?]) could in principle be corrected for 
as well using the value of a, but given the difficulty of achieving 
high accuracy in its determination, that is not recommended. If 
the intensity normalization of the observed spectra is later done 
by normalizing an average profile after the correction to a the- 
oretical or observed reference profile, the loss term will also be 
corrected for automatically by the normalization coefficient. The 
first-order correction is easy to apply and does not introduce nu- 
merical problems when K is given, despite some eventual inter- 
polation in the construction of a discrete 2D kernel. 



2.2. Stray-light correction and spatial deconvolution using a 
known PSF 

To apply the stray-light correction, different approaches are pos- 
sible. One straightforward method is based on the simplified 
Eq. (O. The knowledge of the instrumental PSFinstr also allows, 
however, more detailed methods to be used that take the spa- 
tial dependence of stray light - as given by the second term on 
the right-hand side of Eq. Q - fully into account. For such ap- 
proaches, we describe a first-order and a full correction by a de- 
convolution in the following. 



2.2.2. Fourier deconvolution 

In the Fourier domain, Eq. ([B transforms to 

/obs = K • /true , (9) 

where ~ denotes the Fourier transform and only the static instru- 
mental PSF descri bed by the kernel ^ is considered. It can be 
easily shown (e.g. jDeForest et aDl2009l) that for a known K the 
reverse direction is given by 

/true = • /obs (10) 
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where is the element- wise inverse operator for the Fourier 
transform of K. 

The application of Eq. ([TOl) exceeds the first-order correc- 
tion by taking both the gain and the loss term into account at the 
same time. However, even if the derivation of and the subse- 
quent deconvolution is straightforward, its application can have 
unwanted side eff'ects, such as amplifying noise. We compared 
the Fourier transforr n of the inverse PSF fo r POLIS with the cor- 
responding Fig. 1 of iDeForest et al.l (l2009l) . The resulting curve 
directly resembled the kernel they obtained after regularization 
for the noise suppression. For the present study we therefore as- 
sumed that we do not need to additionally modify the PSF or 
apply a high-frequency noise filtering before the deconvolution. 

3. Stray-light measurements 

The following sections describe the measurements for obtaining 
generic estimates of the stray-light level and of the instrumental 
PSF. We also suggest the best approach to applying the parasitic 
light correction. 

3.1. Parasitic light 

We start with the parasitic light coefficient JS (Eq. O because, 
unlike a, it can be determined directly. From the discussion of 
the source of the parasitic light in the previous section, it fol- 
lows that if the direct illumination of the CCD with sunlight is 
blocked, only the randomly scattered light will be left over. This 
is thus the method of choice in determining J3. We therefore com- 
pared measurements of a spectrum with a fully open aperture 
and a spectrum with the light path blocked close to the camera 
mirror inside of POLIS (Fig. O. The location of the blocking 
prevented the direct illumination of the CCDs with the resolved 
spectrum, but did not aff'ect the stray-light level inside the spec- 
trograph otherwise. Both measurements were corrected for the 
dark current before the analysis. 

To quantify the amount of parasitic light as a fraction of some 
mean intensity of the incoming light, we averaged the intensity 
of the full-aperture spectrum over a pseudo-continuum window 
in the spectrum from 396.383 nm to 396.393 nm as reference, 
and this yielded a mean value of 2716 counts. The average in- 
tensity value of the measurement without direct illumination of 
the CCD was around 148 detector counts. The corresponding im- 
age was homogeneously lit and showed no trace of spectral lines 
which ensures that no dispersed stray light was included in the 
measurement. The ratio of parasitic and reference intensity then 
is about 5 % in the Ca II H channel. The corresponding values for 
the channel at 630 nm were 20380 and 58 counts, respectively, 
yielding a fraction of 0.3 % that is presumably negligible in an 
analysis of the 630 nm data. 

For all following stray-light measurements in the Call H 
channel, we have corrected the spectra for the parasitic light sep- 
arately spectrum by spectrum. We determined the average inten- 
sity in the continuum window in each spectrum, averaging along 
the slit as well, and then subtracted 5% of this average value in 
detector counts from all intensity values of the spectrum. 

For the correction of scientific observations, we suggest us- 
ing a slightly modified approach. From repeated small maps 
taken with about half an hour cadence at the center of the solar 
disk, one can obtain an intensity normalization curve for each 
point of time during the day. Figure [3] shows an example of such 
a curve for 14 August 2009. The first observation of this day 
was a large-area scan on disk center, so its intensity was used 




8.5 9.0 9.5 10.0 10.5 11.0 

UT [h] 



Fig. 3. Intensity normalization curve for the observations of 
14/08/2009. The observed intensity at disk center is given by the 
black crosses, the red line is a polynomial fit to the data points. 



for the calibration curve, too. For each separate scan of a certain 
FOV, one then calculates the average profile over all scan steps, 
first dividing each spectrum by the corresponding value of the 
normalization curve. This removes the temporal trend during the 
scan and provides an average spectrum /locai with an intensity 
only proportional to that of the observed FOV. The parasitic cor- 
rection is then given by the 5%-fraction of the intensity in the 
continuum window of the average profile. To include the change 
in light level during the observation again - because < /locai > at 
a given moment of time scales with it - the normalization curve 
from the start until the end of the scan is normalized to its mean 
value inside that time span. This provides the relative variation 
in the light level, hence of the parasitic correction for a spectrum 
taken at a given moment inside the FOV that entered into the 
calculation of /locai- Multiplying the parasitic correction with the 
relative temporal variation in intensity during the scan yields a 
correction for the parasitic stray light that accounts for both the 
variation in /locai with the telescope pointing and the temporal 
variation in the light level during the day. 

3.2. Spider mounting 

The guiding telescope of the VTT is fed by a small pick-up mir- 
ror located close to the entrance window of the evacuated tele- 
scope tube. The pick-up mirror is fixed to a spider mounting and 
produces a central obscuration in any pupil plane at the VTT. 
We set up an imaging channel to determine the light level inside 
the central obscuration. We deflected the light with the 45-deg 
mirror to an optical rail where we mounted a lens for creating 
a pupil image, a small-band interference filter at 557.6 nm, and 
some neutral density (ND) filters to reduce the intensity to a suit- 
able level for the PCO Sensicam camera used. With this setup, 
we took images of the spider mounting in the pupil plane for 
a CLV run similar to the scientific observations, to measure the 
stray-light level and at the same time to identify a possible varia- 
tion with the telescope pointing. The upper panel of Fig. [4] shows 
the resulting pupil images when moving across the solar disk in 
steps of 100''. We took cuts through the images to determine 
the residual intensity in the shaded areas and normalized each of 
the cuts to its maximum intensity value {lower panel). The pupil 
images show a change in the global light level with the tele- 
scope pointing, but the intensity curves after the normalization 
are nearly identical. The residual intensity in the central obscu- 
ration is always about 4%. The light level at the border of the 
FOV ("tail"), which can only be caused by stray light, is below 
2%. These numbers now, however, refer to a pupil plane, not a 
focal plane. After the reflection on the 45-deg mirror, the light 
passed through only three optical elements (lens, interference 
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Fig. 4. Images of the pupil plane taken with the PCO setup in 
a CLV run at 557 nm. Top panel, pupil images, all shown in an 
identical display range. The offset of the telescope pointing from 
the disk center is given in each subpanel in arcsec. Bottom: in- 
tensities of each pupil image along the cut marked by the red 
vertical line in the image at 



filter, ND filter). Because the spider mounting is located behind 
the entrance window, the stray light caused by the two coelostat 
mirrors and the entrance window is absent. 



3.3. Methods without explicit measurements 

There are three more methods that provide an estimate of stray 
light without requiring explicit measurements. The first is to 
m atch average observed spect ra to an atlas profile, as suggested 
bv lAUende Prieto et al.1 f2004 and lCabrera Solana~era l. (2007). 
A comparison with the presumably stray-light-free spectra of the 
FTS spectral atlas (iKurucz et al.|[l984l : iBrault & Neckellll987h 
yielded the best match for a value of = 20 % for both channels 
of POLIS. The approach is, however, not fully consistent with 
Eq. (O because it implies that all stray light is of parasitic na- 
ture. Instead, the direct measurement of ft of Sect. 13. II excludes 
the high value retrieved with this approach. The approach is also 
aimed at determining the spectral point spread function of the 
optics behind the slit, because it uses an average observed pro- 
file without spatial resolution and the comparison of its shape to 
a reference profile. It therefore measures all spectral resolution 
eff'ects of the spectrograph in addition to the parasitic contribu- 
tion. 

A second method is to compare profiles from the umbra of 
sunspots with only a little intrinsic radiation to an average quiet 
Sun spectrum. iRezaei et~an (l2007l) determined a level of at max- 
imum 12 % of stray light at 396 nm for the sum of a -\- J3 in the 
umbra, which can be taken as the minimum amount of stray light 
in the brighter quiet Sun surroundings. A simultaneous inversion 
of 630 nm and 1565 nm spectra with the SIR code yielded a value 
ofa = 10% of stray fight in the umbra (lBeckll2006ir2008h . 

The last indirect method is the presence of residual polariza- 
tion signal in the telluric line blends of the 630 nm channel. In 
the data reduction of long-integrated (> 20 sees) spectra of the 
630 nm channel, the telluric lines at 630.20 nm and 630.27 nm 
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Fig. 5. POLIS Ca II H spectra with partly illuminated FOV, loga- 
rithmic display range. Left: FOV blocked in F2 (hairline mount- 
ing). Middle: FOV blocked in Fl (filter wheel). Right: slit across 
the solar limb. 



sometimes showed a non-zero polarization signal even if the 
nearby continuum wavelengths showed none. This puzzling be- 
havior could be traced back to the / QUV cross-talk cor- 
rection in combination with the real-time demodulation of the 
Stokes vector in POLIS. This demodulation uses pairs of image 
differences for Stokes QUV, but the addition of all images for 
Stokes /. Parasitic light thus adds up in / but cancels in QUV. 
For a correction of the cross-talk, one then has to use a reduced 
intensity /tme = 4bs - const, in the calculation of the cross-talk 
coefficient cj^quv = G^^Kbs/^tme at continuum wavelengths. 
Subsequently, also cj^quv ' hmc has to be subtracted from QUV. 
This additional correction was used in the data reduction of 
POLIS data in 2003 with const. = 13 = 4 % of 7^, but only for 
the long-integrated data. For shorter integration times, the effect 
was below the noise level and could not be detected. Because the 
stray-light covers in POLIS were modified afterwards as a conse- 
quence of this finding, the number can only serve as an order of 
magnitude estimate. We point out that this effect produces spu- 
rious polarization signals inside spectral lines caused by cross- 
talk from Stokes / even when close-by continuum wavelengths 
are forced to zero polarization. 

3.4. FOV partially blocked in the focal plane (spatial PSF) 

To derive the spatial PSF of POLIS, we took spectra with POLIS 
where part of the FOV in the focal plane was blocked by a metal 
plate. In addition to blocking part of the FOV in the focal planes 
Fl (telescope focus) and F2 (entrance of POLIS), we also took 
spectra with the slit placed across the solar limb. Figure [5] shows 
example spectra of the three cases. The blocking edge is sharper 
for the blocking in F2 inside of POLIS. This is caused by fewer 
optical elements being located behind F2, by the impossibility 
of using the adaptive optics while blocking in Fl and by the fact 
that the metal plate had to be placed at some distance to Fl (^ 1 
cm). All spectra are displayed on a logarithmic scale, otherwise 
the stray light could not be seen at all directly in the images. We 
then averaged the intensity in the same continuum window near 
396.4 nm as before to obtain the spatial intensity variation across 
the blocking edge and the limb, respectively (Fig.O. 

To determine the spatial PSF, we then defined an artificial 
sharp edge located appropriately along the slit, and constructed a 
model for the spatial PSF from a combination of a Gaussian and 
a Lorentzian curve. We convolved the sharp-edge function with 
the kernel and modified the free parameters of the two functions 
used in its generation by trial-and-error until the convolved edge 
function (red line) matched the observed intensity trend across 
the edge for the blocking in Fl. The finally used kernel is shown 
in the right part of Fig.[S]at about x = it drops to basically 

zero at about 10'' distance from its center. The area of the kernel 
is normalized to unity, and its maximum value then is around 
20%. For better visibility, it is shown multiplied by four in the 
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Table 1. Stray-light levels for the VTT/POLIS in various studies. 



Fig. 6. Observed intensity variation along the slit for FOV 
blocked in Fl (black), blocked in F2 (green), and for the slit 
across the limb (blue). The purple lines diix,y = (C, 100%) de- 
note a simulated sharp edge that after convolution with the ker- 
nel at the right yields the red curve. The inlet on the left shows 
a magnification of the stray-light tail region. The derivative of 
the observed intensity variation for the FOV blocked in Fl is 
overplotted over the kernel as an orange line. 



figure. The shape of the kernel can also be cross-checked against 
the observed intensity variation by another method than the con- 
volution of the sharp-edge function. In case of an axisymmetric 
2D kernel K(x,y), it can be shown that the ID kernel K(x) is 
given by 

K{x) = ^{x), (11) 
dx 

when the true intensity corresponds to a step function and 
Ioh^( x) is the observed intensity variation across the step func- 
tion (ICollados & VazquezinWAlBonet et al.lll995h . 

We overplotted the derivative of the observed intensity vari- 
ation for the blocking in Fl over the kernel in the right part of 
Fig.[6l The match is reasonably close, with a slightly lower cen- 
tral value and a broader central core in the derivative. The trial- 
and-error method for defining the kernel thus yielded a good 
match to the theoretical expectation. 

Since the chosen kernel exhibits a sharp drop, most of the lo- 
cal stray light should come from the close vicinity (d < 2") of a 
given point (x, y). The inlet in the left-hand part of Fig. [6] shows 
that, actually, the tail level of the observations is not matched 
well, because the curve of the convolved sharp edge reaches zero 
several arc seconds before the observations. The residual inten- 
sity in the far tail (> 20'') is at maximum 2%. We point out that 
the observations of the solar limb and with the FOV blocked in 
Fl only diff'er slightly for x>^" , which implies that most of the 
stray light is caused by the instrumental PSF and not the tele- 
scope or atmospheric scattering. 

For the application of Eqs. dS]) and ([TOl) , a 2D kernel is 
needed. We used the assumption of axial symmetry to create the 
2D kernel K from the ID version shown in Fig. [Si and divided 
K(x, j) by 2 TT r, with r = -h j^, to maintain the normaliza- 
tion of the area of the ID kernel. 



3.5. Summary of stray-light measurements 

We used diff'erent methods to estimate the amount of stray light 
in POLIS. The upper two sections of Table [T] list the stray-light 
values determined in this work. In the lower section, we added 
the values that have been presented in previous studies for any of 
the two channels of POLIS. The definition of the stray-light frac- 



Measurement 


Wavelength 


Stray-light level 


Pupil plane (PCO) 


557 nm 


4 % (~ 1-2 % tail) 


POLIS (this work) 


Fl blocked 


396 nm 


~ 1-2 % (tail) 


F2 blocked 


396 nm 


~ 1-2 % (tail) 


Limb 


396 nm 


~ 1-2 % (tail) 


Blocked behind grating 


396 nm 


~5 % 


Blocked behind grating 


630 nm 


«; 1 % 


Convolution of FTS 


630/396 nm 


20% 


POLIS 


(other studies) 


Umbral profile ^ 


396 nm 


~ 12 % 


Room light ^ 


396 nm 


< 1 % 


Convolution of FTS ^ 


630 nm 


15 % 


Umbral profile ^ 


630 nm 


~ 10% 


Parasitic light correction 


630 nm 


-4% 



i Rezaei et all (120071) ICabrera Solana et all ([2007h Ib^ 
(2008^ ^"^^ Keyword in POLIS data reduction software, used for 
long-integrated data in 2003. 



tion diff'ers in most cases, but the umbral profiles provide a lower 
limit of 10- 12% for the total stray-light contamination in quiet 
Sun areas without strong spatial intensity gradients. The para- 
sitic light fraction could be accurately determined as 5 % (0.3 %) 
of the intensity for the channel at 396 nm (630 nm) for the data 
taken in 2009, but it presumably was higher in previous years 
before the modification of the stray-light covers inside the spec- 
trograph. The stray-light level at distances above 20'' from any 
lit area is on the order of 2% only, and the spatial PSF is strongly 
peaked with significant contributions only up to about 2''. There 
were no indications that the relative stray light in the telescope 
changes somehow with the telescope pointing in a CLV run. The 
tail of the observed intensity variation across a sharp edge can 
be used as a proxy for the reduction of the stray-light level in 
off-limb spectra. 



4. Application to data 

As a first application to observed data, we compare the perfor- 
mance of the stray-light corrections with and without explic- 
itly using the instrumental PSF on a data set taken near the so- 
lar limb. The application to the limb data set actually yielded 
improved values for the stray-light contribution a because the 
residual intensity beyond the limb provides a good reference for 
the quality of the correction. In the second and third examples, 
we instead test the eff'ects of the deconvolution as a method to 
improve the spatial resolution of observations on spectroscopic 
and spectropolarimetric data, respectively. 

4.1. Example A: application to a limb data set 

To test both the stray-light correction and the spatial deconvo- 
lution, we used a data set that was taken near the solar limb on 
25 August 2009, with an FOV that extended beyond the limb 
(see Fig. [8] below). This off'ers a good opportunity to determine 
the quality of the corrections because the limb is similar to a 
sharp edge and beyond the limb the (pseudo-)continuum inten- 
sity should be zero. The map was taken with 134 steps of Off 3 
step width and an in tegration time of 13 se conds per scan step 
(for more details see iBeck & Rezaeill201 ll) . We applied a vari- 
ation of corrections to the previously flat-fielded data to assess 
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Fig. 7. Comparison of a large-scale average profile multiplied 
by a (thick black) and individual local stray-light profiles calcu- 
lated around a given pixel using the PSF. The labels denote the 
distance of the pixel to the limb in arcsec (negative/positive = 
on/off'-disk). 



the eff'ects of each step on the data. We used the following com- 
binations of corrections: 

1. Parasitic light only (J3). 

2. Parasitic light (J3) with subsequent first-order correction 
(EqM- 

3. Parasitic light (J3) with subsequent Fourier deconvolution 
(Eq.[TO|. 

4. Parasitic light (J3) and stray light (a) as given by Eq. (O 
(more details below). 

Some of the approaches had to be modified slightly for treating 
off'-limb spectra. Method no. 1 serves as a reference of basically 
uncorrected spectra. For the first-order correction in Method 
no. 2, we restricted the sum of Eq. ^ to (x\y') = (x, y) ± (5, 40) 
pixels to avoid using scan steps that were taken at a significantly 
diff'erent time. The Fourier deconvolution in the 3^^ method was 
done according to Eq. ([TOl) by multiplying the Fourier transform 
of the 2D maps of the intensity at each wavelength, 1(A), with 
the inverse of the kernel. 

For Method no. 4, we used a constant value of a = 10% 
as the minimum level of stray light for all spectra on the disk 
and an average local profile /locai calculated from the part of the 
FOV that was most remote from the limb. The intensity beyond 
the limb drops fast, however, with the limb distance (Fig. [6]), so 
subtraction of a fixed a/iocai yielded negative intensity values 
beyond the limb. To avoid this overcompensation, we scaled a 
down with increasing limb distance using the variation in inten- 
sity for the observation with the FOV blocked in Fl (Fig.O. We 
took the observed intensity variation beyond the location of the 
edge (x < O'O and normalized the first point to 10 % to have a 
smooth connection to the correction on the disk. The limb loca- 
tion was set to where the intensity dropped below 5 %. 

It turned out that a good off'-limb correction was only achiev- 
able when the intensity of /locai was increased by a factor of two. 
Because the correction uses a/iocab the factor of two can be at- 
tributed either to a or /locai- Using an average disk center pro- 
file with its intrinsically higher intensity and keeping a at 10 % 
yielded a worse result, because the spectral line blends in the 
wing were located at slightly diff'erent wavelengths than in the 
observed FOV, and it overcompensated for the stray light present 
as well. We therefore used an on-disk correction corresponding 
to a /local with a = 20%, with the off'-limb correction normalized 
to 20 % on the first point. 
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Fig. 8. 2D maps and spectra after the correction. Lower panel: 
OW map in logarithmic display (bottom row), H2R map (middle 
row), spectra (top row) along the cuts marked by white vertical 
lines in the H2R map. Upper panel: magnification of the areas 
marked by black squares in the H2R map (on-disk and at the 
limb). The white rectangles denote the off'-limb area used to de- 
rive the profiles of Fig. [TOl 



The values used in the correction can actually be cross- 
checked with the PSF. Figure |7] shows a /locai used in Method 
no. 4 on the disk, together with some profiles corresponding to 
the gain term of Eq. dS]), calculated locally from the surround- 
ings of each pixel in the FOV weighted with the PSF. The stray- 
light contributions predicted by the gain term match the generic 
^locai well, as long as the pixel, whose surroundings were used, 
was more remote from the limb than 10''. The similarity of the 
stray-light profiles calculated using the PSF to each other and to 
the average profile supports using a single average profile /locai 
in Eq. (|4]) ignoring the PSF, but with a doubled a of 20 %. The 
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Fig. 9. Intensities across the limb in the OW wavelength range 
for an individual spectrum. Solid black line: Method no. 1. 
Triple-dot-dashed orange line: Method no. 2. Dash-dotted blue 
line: Method no. 3. Dotted red line: Method no. 4. 
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Fig. 10. Average off-limb spectra from the area marked by a 
white rectangle in Fig.[8j Line styles and colors are as in Fig.|9j 



profile at 50'' limb distance has a lower intensity because it is 
at the lower border of the FOV, where the full area of the kernel 
could not be used. 

Figure [8] shows the result of applying all four methods (each 
shown in a column of the figure, in the same order as presented 
in the text) to the data set in 2D maps of the FOV and one sample 
spectrum. We derived one map i n a pseudo-contin uum window 
in the outer wing (OW, see, e.g. jBeck eTani2008h of the Call 
H line and one map averaged over the wavelength region of the 
^2R emission peak. The amount of residual off'-limb stray light 
can be clearly seen in the OW map (bottom row) and the spectra 
(third row). Only Method no. 4 removes the off'-limb stray light 
satisfactorily. 

The residual stray-light level beyond the limb is quantified 
in Fig. |9j which shows the intensity variation along the slit in 
the OW map for a randomly chosen scan step. The on-disk light 
level is globally reduced for Methods no. 2 and 4, whereas the 
Fourier deconvolution (no. 3) keeps the same level as the uncor- 
rected spectra (no. 1). This reflects that only with the Fourier 
deconvolution is the loss term of Eq. dS]) taken into account, 
whereas all other methods only correct the gain term. Beyond 
the limb, the stray-light level stays at the value of the uncorrected 
spectra of about 1-2% for the Fourier deconvolution (no. 3), 
whereas the first-order correction (no. 2) and the ad-hoc stray- 
light correction (no. 4) reduced the off'-limb light level. The last 
method yields a correction an order of magnitude better than all 
others, reaching a level of about 0.2 % of residual intensity. 



Fig. 11. Fourier power spectrum of the OW map vs. spatial fre- 
quency. Line styles and colors are as in Fig.[9l 

Table 2. Spatial rms contrast and spectral rms noise in % of 7^. 



method 


1 


2 


3 


4 


contrast 


5.2 


5.8 


6.7 


6.5 


noise 


1.15 


1.18 


1.34 


1.17 



The average profiles shown in Fig. [TOl calculated for the 
off'-limb area marked in the OW map of Fig. [8l also indicate 
the quality of the correction. These profiles allow instant dif- 
ferentiation of under/overcorrected stray light by the shape of 
the line blends in the wing: absorption profiles indicate resid- 
ual stray light, whereas emission profiles for the photospheric 
blends at these heights above the limb indicate an overcorrec- 
tion. Only Method no. 4 yields fully negligible residuals of the 
line blends. For the three other methods, absorption lines are 
seen in the wing. For the emission in the Call H line core - 
basically, the only scientifically interesting quantity beyond the 
limb - Methods no. 2 and 4 yield a similar amplitude, whereas 
the Fourier deconvolution (no. 3) yields nearly the same inten- 
sity as the uncorrected spectra (no. 1). 

Even if Method no. 4 gives the best results for the off'-limb 
stray-light correction, this does not signify that the principle of 
the other methods is wrong. Methods no. 2 and 3 follow the ini- 
tial way in which stray light is created closer than Method no. 4. 
Both approaches, however, only correct the static contribution to 
the stray light caused by the instrumental PSF, within the limits 
of accuracy in its derivation. Their worse performance beyond 
the solar limb implies that in this case both the residual dynam- 
ical part of the PSF beyond the AO correction and the PSF of 
the telescope and the optics upfront of Fl play a role. Because 
the true solar off'-limb intensity outside of the very emission core 
is zero, any residual stray light not accounted for by the instru- 
mental PSF is seen prominently in Figs. [8] or [TOl The ad-hoc 
method (no. 4) instead seems to correct well for both the stray 
light caused by the instrumental PSF, as well as for the addi- 
tional stray-light contributions, but only with a static correction 
without taking fluctuations of the atmospheric scattering into ac- 
count. 

As far as the spatial resolution is concerned, both the first- 
order correction (no. 2) and the Fourier deconvolution (no. 3) 
have led to an improvement, even if it is not very noticeable 
(second and third columns in the top two rows of Fig. [8]). The 
Fourier deconvolution had a stronger eff'ect on the visual im- 
pression. This is verified by the plot of the spatial Fourier power 
for the OW map shown in Fig. [TT] Only the lower part of the 
OW map up to J = 40'' was used. The power at higher spa- 
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Fig. 12. Original {top row) and deconvolved {bot- 
tom row) continuum intensity maps at 777 nm 
of an active region with several pores. Left: full 
FOV. Right: magnification of the region marked 
by a black rectangle. 




tial frequencies is slightly enhanced after the application of the 
first-order correction and significantly enhanced after the Fourier 
deconvolution. The improved spatial resolution is reflected also 
in the root-mean- square (rms) contrast of the OW maps (Table 
|2]). This number, which measures the relative intensity variation, 
has, however, to be taken with caution because any correction 
that reduces globally the intensity automatically increases the 
rms contrast. We also determined the rms variation inside wave- 
length windows of about 10 pm extent in the OW spectral region 
and found it to stay basically at about 1.2% of 7^, with an av- 
erage value of the intensity itself of about 20 % at these wave- 
lengths. The Fourier deconvolution (no. 3) has slightly increased 
the noise. 

4.2. Example B: application to pore observations at 777 nm 

In another observing campaign in November 2009, we took 
identical measurements with a partly blocked FOV with the main 
spectrograph of the VTT in a wavelength range around the O i 
triplet at 777 nm. The corresponding convolution kernel derived 
from the measurement is shown in Appendix|B] The spatial sam- 
pling along the O'.'IS wide slit was O'.'IS in this case. 

We then applied the spatial deconvolution to an observation 
of an active region with several pores, using Eq. ([TOl) without 
any additional stray-light correction. The spatial scanning in the 
observation was done with 201 steps of Qf! 36 width, undersam- 
pling the slit width. Figure [12] shows the full FOV before {upper 
left panel) and after the deconvolution {lower left panel). The 



Fig. 13. 2D maps of the sunspot observations 
at 630 nm. Top/bottom row: original/deconvolved 
spectra. Left to right: Stokes IQUV. The white 
dashed line marks the location of the cut shown 
in Fig. [T4j the crosses and the corresponding 
numbers label the locations of the profiles shown 
in Fig. [15] The horizontal black/white lines at the 
upper and lower border of the FOV are caused by 
the hairlines of POLIS. 

seeing was moderate to bad during the scan, producing frequent 
failures of the AO system, with jumps of the FOV as well. The 
magnification of two of the pores {right panels), however, shows 
that the deconvolution has not only increased the contrast by a 
reduction of the intensity, but also leads to an improvement in 
the sharpness of the structures, e.g., for the boundaries of the 
darkest parts of the pores. 

4.3. Example C: application to polari metric observations 

To test the eflTect of the deconvolution by Eq. ([TOb on polarimetric 
data, we applied it to observations taken with POLIS in its red 
channel at 630 nm. Assuming that the PSF derived for the blue 
channel is characteristic of the instrument, we did not determine 
a new PSF but used the same one as before. Appendix |B] shows 
that applying this kernel to a knife-edge observation at 630 nm 
also provides a satisfactory match. The observation we discuss 
here consists of a scan with 150 scan steps of 0'.'3 step width 
across an active region (Fig. [13]) done on 6 July 2005, with AO 
correction. The integration time was ten seconds per scan step. 
In this case, the Fourier deconvolution was applied separately to 
each wavelength and Stokes parameter. 

In the visual impression of Fig. \T3\ the improvement in the 
spatial resolution is prominent, both in the continuum intensity 
map {left) and the maps of the absolute wavelength-integrated 
polarization signal, J \QUV\/IdA. The umbral dots inside the 
umbra of the bigger spot can be seen well in the deconvolved 
data set and the filamentary structure of the penumbra of the 
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20 30 
arcsec 



Fig. 14. Cut through the intensity maps of Fig. [131 Thick black: 
original data. Thin grey, deconvolved data. The upper panel 
shows a magnification of the central region where the cut crossed 
a light bridge. 




630.2 630.3 630.4 
A [nm] 

Fig. 15. Stokes V spectra from the sunspot observation. Red: de- 
convolved spectra. Black: original spectra. The numbers in the 
upper left of each panel refer to the labels in Fig. [13] that indi- 
cate the locations of the spectra in the 2D maps of the sunspot 
observation. 



smaller spot is much clearer. The rms contrast in the continuum 
image increased from 3.7% to 5.3%. Figure [T4l demonstrates 
that the increase in the rms contrast is not only caused by the 
removal of the stray-light contamination. It shows a cut through 
the umbra of the larger sunspot where a light bridge was located. 
The full width at half maximum of the light bridge is slightly re- 
duced and the intensity values to the left {x ~ 26'0 and right 
{x ~ 30.5'0 of the central region decrease to the benefit of a 
more concentrated distribution with a stronger intensity peak in 
the deconvolved data, compared to the original data. 

The polarization signal was also usually amplified by the 
deconvolution. Figure [T5] shows four Stokes V profiles before 
and after the deconvolution; taken from the locations marked in 
Fig. [131 On average, the Stokes V amplitudes increased relatively 
by about 5 % of their previous value. These spectra, however, 
also clearly show the drawback of the Fourier deconvolution: the 




a [%] 



Fig. 16. Histograms of the local value of a in a QS map {thick 
black), a sunspot map (red), and the limb data set (thin blue). 



noise level is amplified as well. In Stokes /, the rms noise in a 
continuum window increased from 0.05 % to 0.07 % at an aver- 
age value of about 15. 000 counts. In Stokes V, where the noise 
level is more critical in the analysis of the data, the rms noise was 
0.05 % of Ic before and 0.12 % of after the deconvolution, re- 
spectively. The rms noise in Stokes V is still acceptabl e, even af- 
ter the increase by a factor of 2.4. For comparison, Steiner et al.l 
(2010) found an increase of the rms noise from 0.1 % to 0.2- 
0.3 % in the deconvolution of IMaX data. 



5. The relation between PSF and (local) stray light 

In the stray-light correction of the off'-limb spectra, it turned out 
by trial-and-error that a value of a = 10% is definitely too 
low. The finally used stray-light level that yielded a good cor- 
rection corresponds to using a = 20%. The integration of the 
local stray-light contribution using the PSF finally returned the 
same level (Fig. O. As seen in Eq. (j?)), the PSF and the (lo- 
cal) stray light are mutually related to each other in a way that 
makes a separation between them difficult. To cross-check that 
the approach of using a single profile averaged over the full FOV 
instead of the explicit local version given by the PSF is reason- 
able, we calculated the local stray light using the gain term of 
Eq. © for all pixels of three data sets: the limb (@396nm) and 
sunspot (@630nm) observations used in the previous section, 
and a large-area map taken in quiet Sun (QS) on disk center 
(@396 nm, not shown here). We then determined the intensity of 
the local stray-light profiles relative to the average profile of the 
full FOV because the number can be used as first-order prox}Qof 
the true a. Figure [16] shows the histograms of the relative fre- 
quency of occurrence of a given value of a. In QS regions, a is 
about 24 %. For the limb and the sunspot data sets it is also gen- 
erally at or above 20%; however, it can in some instances fall 
below 15 % (see the tail of the corresponding distributions in the 
figure). Throughout the umbra of the large sunspot of Fig. \T3\ a 
was about 8-10%. 



6. Discussion 

The various stray-light measurements are basically consistent 
with the formulation of the stray-light problem in Eqs. (0, ©, 
and (O. For the special case of POLIS, the stray light contains 
a significant contribution of spectrally undispersed light created 



1 In Eq. HI), /obs is used because /tn 
slightly overestimated. 



is unknown; a is therefore 
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inside the instrument itself that could be absent in other spec- 
trographs. For the spectrally resolved stray light that is created 
by scattering in the light path upfront of the grating, we found 
a lower limit of about 10 % from umbral profiles and a value of 
about 20 % in the quiet Sun, consistently from a limb data set 
and an application of the PSF to calculate the local stray light. A 
deconvolution with the instrumental PSF alone, however, does 
not provide a good correction for off'-limb spectra, presumably 
because it only covers the optics behind the telescope focus and 
includes neither the atmospheric scattering nor the telescope PSF 
of the optics upstream of the focal plane. For off'-limb data, addi- 
tional corrections are needed that have eventually to be adjusted 
to the specific data set used, but a residual stray-light level be- 
low 0.5 % can be achieved. For our data, we used the observed 
intensity variation across a sharp edge with only a minor ad-hoc 
adjustment to satisfactorily remove the off'-limb stray light. 

We determined the PSF from observations with a partly 
blocked FOV, for the two channels of POLIS and one obser- 
vational run with the main spectrograph of the VTT at 777 nm. 
For the latter and the 630 nm channel of POLIS, the tail region 
of the stray light far away from the blocking edge was only par- 
tially reproduced, i.e., the derived kernel yields less stray light 
than actually observed. This is possibly related to our basically 
ID approach using individual slit spectra. A spatial scan of the 
blocking edge and the determination of a 2D PSF may improve 
the result in the tail region because in the ID case only the con- 
tamination along the slit is included, and not the "lateral" stray 
light. Determination of the instrumental PSF from observations 
with a half-blocked FOV is not limited to slit- spectrograph data 
but can equally be done for all kind of ID or 2D instruments. 
It partly allows one to avoid hav i ng to u se purely a d -hoc c on- 
structed PSFs as in iPereira et all (l2009h . iJoshi et al.l (l201lh . or 
iRutten et al.l (1201 ll) . The derivation of the PSF from explicit ob- 
servations is more solid than the indirect method of a comparison 
between observed an d synthetic spectra from simulations (e.g., 
IScharmer et al.ll20TTl) , where spatial and spectral resolution ef- 
fects (cf. Sect. 13.31) can get mixed up. 

The instrumental PSF can be used to deconvolve the ob- 
served spectra post-facto for both a correction of the stray light 
and an improvement of the spatial resolution. We tested the de- 
convolution on three diff'erent data sets, two spectroscopic and 
one spectropolarimetric observation. From the visual impres- 
sion, the spatial resolution has improved in all cases beyond a 
simple increase in contrast owing to the global reduction of in- 
tensity because of the subtraction of the stray-light contribution. 
The noise level was increased by the deconvolution by up to a 
factor of about 2, while our PSF seems to be well behaved in gen- 
eral with respect to the noise amplification, presumably because 
it is sampling-limited in all cases and does not cover the high 
spatial frequencies down to the theoretically achievable diff'rac- 
tion limit where the noise contribution to the Fourier power also 
becomes important. For the polarimetric observation, we point 
out that we used the instrumental PSF to deconvolve an obser- 
vation taken five years before the PSF measurement. Even if this 
still led to a clear improvement, the PSF measurement should 
be done closer in time. Before such a deconvolved data set can, 
however, be used for a scientific investigation, the eff'ects of the 
deconvolution on the physically relevant quantities like the net 
circular polarization or the LOS vel ocities would have to be in - 
vestigated in more detail (see, e.g.. iPuschmann & Beckll2QT]]) . 
This should include tests with synthetic observations that are 
convolved with the known PSF. 

There are, however, more reasons that the deconvolution 
works well with the data sets used. The PSF is sharply peaked. 



therefore pixels with a separation of more than about 2-3'' from 
a spatial location do not enter strongly. This automatically re- 
stricts the correction to the close surroundings of a given pixel, 
which is necessary for slit- spectrograph data with their sequen- 
tial spatial scanning over a permanently evolving solar scene. A 
second reason is that all observations were obtained using AO 
correction, therefore the time- variant part of the instantaneous 
PSF is already largely compensated for. The third reason is that 
all observations initially have a high signal-to-noise (S/N) ratio 
because of the comparably long integration times of ten seconds 
or more. Finally, the observations were also spatially oversam- 
pled for the data used here; i.e., the spatial sampling along and 
perpendicular to the slit was below the eff'ective resolution by a 
factor of two or more. 

7. Conclusions 

Observations with a partly blocked FOV in the first accessible 
focal plane allow one to estimate the (static) instrumental PSF 
downstream of the focal plane, in case no observations of a plan- 
etary transit are available, but it seems generally to be recom- 
mended to use a 2D FOV to determine the PSF. The PSF ob- 
tained allows one to determine the local stray-light contamina- 
tion of a pixel from its close surroundings. For a homogeneously 
lit area - as is common in the quiet Sun - the explicit calculation 
by the PSF can be exchanged by an average stray-light contri- 
bution proportional to the average profile of the full FOV; in the 
case of sunspot observations, the stray-light correction should be 
calculated explicitly using the PSF because of the large intensity 
gradients. The PSF and the corresponding stray-light level can 
be cross-checked with observations near the solar limb, where 
continuum wavelength windows provide an excellent intensity 
reference. For ground-based observations taken with AO correc- 
tion at the VTT, with a field stop before the first focal plane, the 
instrumental PSF seems to be dominating over the contributions 
from the telescope and the Earth's atmosphere. 

The PSF can also be used for a spatial deconvolution using 
a Fourier method. This seems a valid option for correcting slit- 
spectrograph observations post-facto for the static contributions 
of the instrumental PSF. 
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Appendix A: Derivation of first-order correction 

Solving Eq. d?]) for /tme yields at first the recursive equation 



1 



where r = (x - x\y - /). 



/obsUj)-^/^(r)/trueU',/) 



x',y' 



(A.1) 




Fig.B.l. Derivation of the PSF at 111 nm. Thick black: ob- 
served intensity. Dashed: edge function. Thin red: convolved 
edge. Upper panel: observed intensity across the blocked edge. 
The convolution kernel is shown on the right; for better visibility 
it was multiplied by 2. Lower panel: "tail" region. 



Exchanging the term /tme(-^', with the relation given by 
Eq. (lA.ll) gives 



1 



x',y' x",y" 



(A.2) 



With 



I -\- a for a <^ 1 and neglecting all terms 
oc a,a^, K a, K^, one obtains 

/tmeU y) = /obsU, y)-YjK(x-x\y- y )/obs(x', y ) . (A.3) 

x',y 

Appendix B: Convolution kernels at 777 nm and 
630 nm 

Convolution kernel at 777 nm Figure iBJl shows the observed 
intensity across the location of the blocking edge at 777 nm, to- 
gether with the convolution kernel used to match the convolved 
edge function to the observation (upper panel). The kernel has 
more extended wings than the one of POLIS, but we were still 
not able to perfectly reproduce the tail region for distances over 
about 15'' from the location of the edge (lower panel). Further 
increasing the wing contribution led to a mismatch on the other 
side of the edge in the lit area because the intensity of the con- 
volved edge curve gets reduced more. Up to about W distance 
from the edge, the observations are reproduced well by the con- 
volved edge. 

Convolution kernel at 630 nm For the 630nm channel of 
POLIS, we resampled the kernel determined at 396 nm to the 
two times finer spatial sampling in the red channel by interpo- 
lation. Figure lR2l shows that an application of this kernel to the 
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Fig. B.2. Test of the PSF at 630 nm. Thick black: observed in- 
tensity. Dashed: edge function. Thin red: convolved edge. The 
convolution kernel is shown on the left; it was multiplied by 8 
for better visibility. 

edge function reproduces the intensity in the 630 nm channel for 
the observation with a half-blocked FOV satisfactorily, even if 
in the tail region the match is worse than for 396 nm. Because of 
the resampling, the amplitude of the kernel is halved compared 
to Fig.[6l 
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